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We present a measurement of the top quark mass with the Matrix Element method in the lep- 
ton+jets final state. As the energy scale for calorimeter jets represents the dominant source of 
systematic uncertainty, the Matrix Element likelihood is extended by an additional parameter, 
which is defined as a global multiplicative factor applied to the standard energy scale. The top 

quark mass is obtained from a fit that yields the combined statistical and systematic jot energy 
scale uncertainty. Using a data set of 370 pb~^ taken with the DO experiment at Run II of the 
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Fermilab Tevatron Collider, the mass of the top quark is measured using topological information to 
be: 

mf+j''*"(topo) = 169.2l?:5 (stat. + JES) tH (syst.) GeV , 
and when information about identified b jets is included: 

m^+j"*"(6-tag) = ITO.alts (stat. + JES) H j (syst.) GeV . 

The measurements yield a jet energy scale consistent with the reference scale. 



PACS numbers: 14.65.Ha, 12.15.Ff 

I. INTRODUCTION 

The origin of mass of the elementary fermions of the 
standard model is one of the central questions in par- 
ticle physics. Of the six known quark flavors, the top 
quark is unique in that its mass can be measured to the 
percent-level with the current data of the Fermilab Teva- 
tron Collider. There is particular interest in a precision 
measurement of the top quark mass because of its domi- 
nant contribution in loop corrections to electroweak ob- 
servables such as the p parameter. Within the standard 
model, a precise determination of the top quark mass 
in combination with existing electroweak data can place 
significant constraints on the mass of the Higgs boson [1] . 

To date, proton-antiproton collisions at the Fermilab 
Tevatron Collider provide the only possibility to produce 
top quarks. During Run I of the Tevatron in the 1990s, 
at a proton-antiproton center-of-mass energy of = 
1.8 TcV, the top quark was discovered by the CDF and 
DO [2] experiments, and its mass was measured [3]. Since 
the beginning of Run II in 2002, the Tevatron is running 
with an increased luminosity at a center-of-mass energy 
of ^/s = 1.96 TeV. The CDF experiment has recently 
published a measurement of the top quark mass using 
Run II data [4]. The first DO measurement of the top 
quark mass at Run II is described in this paper. 

Top quarks are produced in proton-antiproton colli- 
sions either in pairs (production of tt pairs via the strong 
interaction) or singly (via the electroweak interaction). 
Only the first process has been observed so far and is used 
to measure the top quark mass. In the standard model, 
the top quark essentially always decays to a 6 quark and 
a real W boson. The topology of a tt event is therefore 
determined by the subsequent W boson decays. The so- 
called lepton-hjets topology, where one W boson decays 
to an electron or muon and the corresponding neutrino 
while the other decays hadronically, allows the most pre- 
cise experimental measurement of the top quark mass. 
These events are characterized by an energetic, isolated 
electron or muon (charge conjugate modes are implicitly 
included throughout this paper), missing transverse en- 
ergy relative to the beamline from the neutrino, and four 
energetic jets. 

This paper describes a measurement of the top quark 
mass with the DO detector, using lepton-l-jets events from 

370 pb"-'^ of data collected during Rim II of the Fermilab 
Tevatron Collider. To make maximal use of kinematic 



information, the events selected are analyzed with the 
Matrix Element method. This method was developed by 
DO for the Run I measurement of the top quark mass [5] 
and led to the single most precise measurement during 
Run I. For each event, a probability is calculated as a 
function of the top quark mass that this event has arisen 
from tt production. A similar probability is computed for 
the main background process, which is the production of 
a leptonically decaying W boson produced in association 
with jets. The detector resolution is taken into account in 
the calculation of these probabilities. The top quark mass 
is then extracted from the joint probability calculated for 
all selected events. To reduce the sensitivity to the energy 
scale of the jets measured in the calorimeter, the Matrix 
Element method has been extended so this scale can be 
determined simultaneously with the top quark mass from 
the same event sample [4, 6]. 

The paper is organized as follows. Section II gives 
a brief overview of the DO detector. The event recon- 
struction, selection, and simulation are discussed in Sec- 
tion III. A detailed description of the Matrix Element 
method is given in Section IV. The top quark mass fit 
is described in Sections V and VI for the analyses be- 
fore and after the use of 6-tagging information, and Sec- 
tion VII lists the systematic uncertainties. Section VIII 
summarizes the results. 



II. THE DO DETECTOR 

We use a right-handed coordinate system whose origin 
is at the center of the detector, with the proton beam 
defining the positive z direction. The DO detector con- 
sists of a magnetic central-tracking system, comprising 
a silicon microstrip tracker (SMT) and a central fiber 
tracker (CFT), both located within a 2 T superconduct- 
ing solenoidal magnet [7]. The SMT has « 800,000 in- 
dividual strips, with typical pitch of 50 — 80 /xm, and a 
design optimized for tracking and vertexing capability at 
pseudorapidities of jryj < 2.5. The system has a six- barrel 
longitudinal structure, each with a set of four layers ar- 
ranged axially around the beam pipe, and interspersed 
with 16 radial disks. The CFT has eight thin coaxial 
barrels, each supporting two doublets of overlapping scin- 
tillating fibers of 0.835 mm diameter, one doublet being 
parallel to the collision axis, and the other alternating by 
±3° relative to the axis. Light signals are transferred via 
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clear fibers to solid-state photon counters (VLPC) that 
have ~ 80% quantum efficiency. 

Central and forward preshower detectors located just 
outside of the superconducting coil (in front of the 
calorimctry) are constructed of several layers of extruded 
triangular scintillator strips that are read out using 
wavelength-shifting fibers and VLPCs. The next layer of 
detection involves three liquid- argon/uranium calorime- 
ters: a central section (CC) covering jryj up to w 1.1, 
and two end calorimeters (EC) that extend coverage to 
Ir^l « 4.2, all housed in separate cryostats [8]. In addition 
to the preshower detectors, scintillators between the CC 
and EC cryostats provide sampling of developing showers 
at 1.1 < \r]\ < 1.4. 

A muon system [9] is located beyond the calorimctry 
and consists of a layer of tracking detectors and scintilla- 
tion trigger counters before 1.8 T iron toroids, followed by 
two similar layers after the toroids. Tracking at \r]\ < 1 
relies on 10 cm wide drift tubes [8], while 1 cm mini-drift 
tubes are used at 1 < |r/| < 2. 

Luminosity is measured using plastic scintillator arrays 
located in front of the EC cryostats, covering 2.7 < \ri\ < 
4.4. Trigger and data acquisition systems are designed to 
accommodate the high luminosities of Run II. Based on 
preliminary information from tracking, calorimetry, and 
muon systems, the output of the first level of the trigger 
is used to limit the rate for accepted events to ~ 2 kHz. 
At the next trigger stage, with more refined information, 
the rate is reduced further to « 1 kHz. These first two 
levels of triggering rely mainly on hardware and firmware. 
The third and final level of the trigger, with access to 
all the event information, uses software algorithms and a 
computing farm, and reduces the output rate to w 50 Hz, 
which is written to tape. 



III. EVENT RECONSTRUCTION, SELECTION, 
AND SIMULATION 

This paper describes the analysis of 370 pb^^ of data 
taken between April 2002 and August 2004. Events con- 
sidered for the analysis must initially pass trigger condi- 
tions requiring the presence of an electron or muon and 
a jet. In the e-|-jets trigger, an electron with transverse 
momentum pt > 15GeV within \r]\ < 1.1 is required. In 
addition, a jet reconstructed using a cone algorithm with 
radius ATZ = ^ (Ar/)^ -I- (Ai/))^ = 0.5 is required, with a 
minimum jet transverse energy, Et-, of 15 or 20GeV de- 
pending on the period of data taking. In the /i+jets trig- 
ger, a muon detected outside the toroidal magnet (corre- 
sponding to an effective minimum momentum of around 
3 GeV) is required, along with a jet with transverse en- 
ergy Et of at least 20 or 25 GeV, depending on the data 
taking period. 

The offline reconstruction and selection of the events 
is described in detail in the following sections. The kine- 
matic selection criteria are also summarized in Table I. 
A total of 86 e-|-jets and 89 /u-|-jets events are selected. 



, J, , St > 20 GeV bl < 1.1 (electrons) 

charged lepton ^oa/^t/ ^on) ^ 

PT>20GeV |r;| <2.0 (muons) 

exactly 4 jets £t > 20 GeV | | < 2 . 5 

missing transverse energy ^ > 20 GeV 

TABLE I: A summary of the kinematic event selection. In 
addition, quality and isolation criteria are applied. 



A. Charged Lepton Selection 

Candidates for a charged lepton from W decay are re- 
quired to have a transverse momentum of at least 
20 GeV and must be within \r]\ < 1.1 for electrons and 
1 77 1 < 2.0 for muons. In addition, charged leptons have to 
pass quality and isolation criteria described below. 

Electrons must deposit at least 90% of their energy in 
the electromagnetic calorimeter within a cone of radius 
ATZ = 0.2 around the shower axis. The transverse and 
longitudinal shower shapes must be consistent with those 
expected for an electron, based on Monte Carlo simula- 
tion, with efficiencies corrected for observed differences 
between data and Monte Carlo. A good spatial match 
of the reconstructed track in the tracking system and 
the shower position in the calorimeter is required. Elec- 
trons must be isolated, i.e., the energy in the hollow cone 
0.2 < ATZ < 0.4 around the shower axis must not exceed 
15% of the electron energy. Finally, a likelihood is formed 
by combining the above variables with information about 
the impact parameter of the matched track relative to 
the primary interaction vertex, the number of tracks in 
a cone with radius ATZ = 0.05 around the electron can- 
didate, the Pt of tracks (excluding the track matched to 
the electron) within a cone with radius ATZ = 0.4, and 
the number of strips in the central preshower detector 
associated with the electron. The value of this likelihood 
is required to be consistent with expectations for high-py 
isolated electrons. 

For each muon, a match of muon track segments in- 
side and outside the toroid is required. The timing infor- 
mation from scintillator hits associated with the muon 
must be consistent with that of a particle produced in 
the pp collision, thereby rejecting cosmic rays. A track 
reconstructed in the tracking system and pointing to the 
event vertex is required to be matched to the track in 
the muon system. The muon must be separated from 
jets, satisfying ATZ{ix, jet) > 0.5 for all jets in the event. 
Finally, the muon must pass an isolation criterion based 
on the energy of calorimeter clusters and tracks around 
the muon: The calorimeter transverse energy in the hol- 
low cone 0.1 < ATZ < 0.4 around the muon direction is 
required to be less than 8% of the muon transverse mo- 
mentum, and the sum of transverse momenta of all other 
tracks within a cone of radius ATZ = 0.5 around the muon 
direction must be smaller than 6% of the muon pT- 
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B. Jet Reconstruction and Selection 

Jets are defined using a cone algoritlini with radius 
A7^ = 0.5. They arc required to have pscudorapidity 
\r]\ < 2.5. Calorimeter cells with negative energy or with 
energy below four times the width of the average elec- 
tronics noise are suppressed (unless they neighbor a cell 
of high positive energy, where the threshold is lowered 
by a factor of two) in order to improve the calorimeter 
performance. In the reconstruction, jets are considered 
only if they have a minimum raw energy of 8 GeV. Jets 
must then pass the following quality requirements: 

• the energy reconstructed in the electromagnetic 
part of the calorimeter must be between 5% and 
95% of the total jet energy; 

• the fraction of energy in the outer hadronic 
calorimeter must be below 40%; 

• the energy ratio of the most and second most en- 
ergetic calorimeter cells in the jet must be below 
10; 

• the most energetic calorimeter cell must not contain 
more than 90% of the jet energy; 

• the jet is required to be confirmed by the indepen- 
dent trigger readout; and 

• jets within ATZ < 0.5 of an isolated electromagnetic 
object (electron or photon) with Et > 15 GeV re- 
constructed in the calorimeter are rejected. The 
electromagnetic objects used here are obtained 
with a selection similar to the electron selection 
described in Section III A, but without the require- 
ment of a track match or a cut on the likelihood. 

The analysis is restricted to events with exactly four 
jets; these four jets must each have Et > 20 GeV af- 
ter jet energy scale correction, which is described below. 
The motivation for the requirement of four jets is that 
for each event, a signal probability Pgig is calculated us- 
ing a leading-order matrix element for fl production, as 
described below in Section IV C. Decays in which addi- 
tional radiation is emitted as well as ti pairs produced in 
association with other jets are not modeled in the prob- 
ability. 



The corrections are applied to account for several effects: 

• Energy Offset Eon'- Energy in the clustered cells 
which is due to noise, the underlying event, multi- 
ple interactions, energy pile-up, and uranium noise 
lead to a global offset of jet energies. This offset 
EoS is determined from energy densities in mini- 
mum bias events. 

• Showering Corrections Ccone: A fraction of the 
jet energy is deposited outside of the finite-size jet 
cone. Jet energy density profiles are analyzed to 
obtain the corresponding correction Ccone- 

• Calorimeter Response i?cai: Jets consist of 
different particles (mostly photons, pions, kaons, 
(anti-)protons, and neutrons), to which the 

calorimeter response differs. Furthermore, the en- 
ergy reponse of the calorimeter is slightly nonlin- 
ear. The response i?cai is determined from 7-l-jets 
events by requiring transverse momentum balance. 
The photon energy scale is assumed to be identical 
to the electron scale and is measured independently 
using Z ^ ee events. 

Note that E^°^'' is not the parton energy: the parton may 
radiate additional quarks or gluons before hadronization, 
which may or may not end up in the jet cone. The rela- 
tion between the jet and parton energies is parameterized 
with a transfer function, see Section IV B. The jet energy 
scale is determined separately for data and Monte Carlo 
jets. The scale depends both on the energy of the jet and 
on the pscudorapidity. All jet energies in data and Monte 
Carlo events are corrected according to the appropriate 
jet energy scale, and these corrections are propagated to 
the missing transverse energy, see Section HID. 

The uncertainty on the jet energy scale was the domi- 
nating systematic uncertainty on most previous measure- 
ments of the top quark mass. To reduce this systematic 
uncertainty, information from the jets arising from the 
hadronic W decay can be used to determine an overall 
jet energy scale factor, JES, simultaneously with the top 
quark mass. A value of JES ^ 1 means that all jet ener- 
gies need to be scaled by a factor of JES relative to the 
reference scale described above. 



D. Missing Transverse Energy 



C. Jet Energy Scale 

The measured energy -Ejl™ of a reconstructed jet is 
given by the sum of energies deposited in the calorime- 
ter cells associated with the jet by the cone algorithm. 
The energy E^°" of the jet before interaction with the 
calorimeter is obtained from the reconstructed jet energy 
as 

preco TP a 

^corr _ -^jet -^off 

-^cal^cone 



Neutrinos can only be identified indirectly by the im- 
balance of the event in the transverse plane. This imbal- 
ance is reconstructed from the vector sum of all calorime- 
ter cells with significant energy (cf. Section IIIB). The 
missing transverse energy is corrected for the energy scale 
of jets and for muons in the event. Only events with 
IpT > 20 GeV are considered. 

In addition, a cut on the difference A(j) between the 
azimuthal angles of the lepton momentum and the miss- 
ing transverse energy vector is imposed to reject events 
in which the transverse energy imbalance originates from 
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a poor Icpton energy measurement. This requirement 
depends on the scalar value of If^T and is 

A</,(e,^)> 0.777-^^ (2) 
in the e+jets channel and 

Act> (m, ^) < tt - 0.277 (^1 - (3) 
in the /U+jets channel. 

E. b Jet Identification 

A tt event contains two 6 jets, while jets produced in 
association with W bosons predominantly originate from 
light quarks or gluons. The signal to background ratio is 
therefore significantly enhanced when requiring that one 
or more of the jets be idciitificxi as b jets (6-tagged). DO 
developed a lifetime based /^-tagging algorithm referred to 
as SVT [10]. The algorithm starts by identifying tracks 
with significant impact parameter with respect to the 
primary vertex. Only tracks that are displaced by more 
than two standard deviations are considered. The algo- 
rithm then requires that these tracks form a secondary 
vertex displaced by more than seven standard deviations 
from the primary vertex. For each track participating 
in secondary vertex reconstruction its impact parameter 
must have a positive projection onto the jet axis. Tracks 
with a negative projection appear to originate from be- 
hind the primary vertex, which is a sign of a mismeasure- 
ment. Tracks with negative impact parameter are used 
to quantify the mistagging probability. 

The performance of the SVT algorithm is extensively 
tested on data. The 6-tagging efficiency is verified on a di- 
jet data sample whose b jet content is enhanced by requir- 
ing that one of the jets be associated with a muon. The 
distribution of the transverse momentum of the muon 
relative to the associated jet axis is used to extract the 
fraction of b jets before and after tagging. The proba- 
bility to tag a light quark jet (mistag rate) is inferred 
from the rate of secondary vertices with negative impact 
parameter, corrected for the contribution of heavy flavor 
jets to such tags and the presence of long-lived particles in 
light quark jets. Both corrections are derived from Monte 
Carlo simulation. Both the 6-jet tagging efficiency ejet(&) 
and the light-jet tagging rate ejet(w, d, s, g) are pa- 
rameterized as functions of the transverse jet energy and 
pseudorapidity. The efficiency ejct(c) to tag a c quark 
jet is estimated based on the Monte Carlo prediction for 
the b to c-jet tagging efficiency ratio. These parameteri- 
zations are used to predict the probability for a jet of a 
certain flavor to be tagged. 



F. Simulation 

Large samples of Monte Carlo simulated events are 
used to determine the detector resolution, to calibrate 
the method, and to cross-check the results for the top 
quark mass obtained in the data. The alpgen [11] event 
generator is used for both signal and background simu- 
lation. The hadronization and fragmentation process is 
simulated using PYTHIA [12]. Signal tt events arc sim- 
ulated for top quark mass values of 160, 170, 175, 180, 
and 190 GeV. The main background is from jets 
events. 

All simulated events are passed through a detailed sim- 
ulation of the detector response and are then subjected 
to the same selection criteria as the data. The probabil- 
ity that a simulated event would have passed the trigger 
conditions is calculated, taking into account the relative 
integrated luminosities for which the various trigger con- 
ditions were in use. This probability is typically between 
0.9 and 1 and is accounted for when simulated events are 
used in this analysis. 

Background from QCD multijet processes has not been 
generated in the simulation; instead, events that pass 
a selection with reversed isolation cuts for the charged 
lepton have been used to model this background. 

G. Sample composition 

Even though the Matrix Element method yields the ti 
content /top of the selected data sample together with 
the top quark mass and jet energy scale, an independent 
estimate of the sample purity is obtained using a topolog- 
ical likelihood discriminant, as described in this section. 
This result for the sample composition does not directly 
enter in the top quark mass fit; it is used only 

• to obtain the relative normalization of the signal 

and background probabilities as described in Sec- 
tion IV D — this allows fitting the sample purity 
without large corrections to the result — and 

• to choose the sample purity in ensemble tests (cf. 
Sections V B, V C, and VI B) according to the sam- 
ple composition in the data — in order to compute 
the expected fit uncertainties. 

The relative contributions of tt, W^+jets, and QCD 
multijet events to the selected data sample are deter- 
mined before b tagging is applied. A likelihood discrimi- 
nant based on topological variables is calculated for every 
selected event. The technique is the same as described 
in [13] . A fit to the observed distribution yields the frac- 
tions of tt, VK+jets, and multijet events in the data sam- 
ple, separately for e+jets and /x-|-jets events. The fits 
are shown in Figure 1, and the results are summarized 
in Table II. Note that because of differences between the 
selection criteria of e-|-jets and /i-|-jets events (most no- 
tably, the 1 77 1 requirement, but also the criteria for select- 
ing isolated leptons), the numbers of selected tt, W+]eis, 
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FIG. 1: Distribution of the topological likelihood for the 
370 pb~^ DO Run II data sample. The distribution for e+jets 
events is shown in plot (a) and for /i+jets events in plot (b). 
The points with error bars indicate the data, and the fit- 
ted fractions of ti events (open area), VK+jets events (di- 
agonally hatched), and QCD multijet events (horizontally 
hatched area) are superimposed. 



and QCD events in the two channels are not expected to 
be equal. 

The relative contributions from background events 
with a W boson and four jets with different flavor com- 
position are estimated using the ALPGEN generator. The 
fractions /$ of each of the six flavor configurations $ = 
jjjj, bbjj, ccjj, {bb)jjj, {cc)jjj, and cjjj are listed in 
Table III [10]. The symbol j denotes a light jet not con- 
taining a charm or bottom quark, and the symbols (bb) 
and (cc) refer to situations where two heavy flavor quarks 
end up in the same jet. These fractions are obtained 
without a 6-tagging requirement; in the 6-tagging anal- 
ysis where the event sample is divided into three classes 
according to the number of 6-tagged jets per event (cf. 
Section VIA), the individual numbers for each of the 



channel 


Ncvts 


top 


/■topo 
J top 




rtopo 

/qcd 


e-|-jets 


86 


40.6 1 l:t 


(47.211°;^) 


% 


(17.6 1 li) % 


/i-l-jets 


89 


25.8 t li 


(29.01 ID 


% 


(5.11 S;^)% 


^-l-jets 


175 


66.4 


(37.91 1:1) 


% 


(11.3 ±1.2)% 



TABLE II: Composition of the e-|-jets, ^+jets, and £+]ets 
data samples, estimated with the topological likelihood tech- 



nique. The fractions are constrained as /, 

/topo -1 
QCD — ^■ 



topo 
top 



topo 
W+jcts 



Contribution W + > 4: jets 



Wbbjj 

Wccjj 

W{bb)jjj 

W{cc)jjj 

Wcjjj 

Wjjjj 



(2.72 ± 0.11)% 

(4.31 ± 0.20) % 

(2.70 ± 0.15)% 

(4.69 ± 0.36) % 

(4.88 ± 0.17)% 

;80.71 ± 0.43) % 



TABLE III: Fractions /$ of different flavor subprocesses con- 
tributing to the W-|-jets sample. 



three separate classes are significantly different. 



IV. THE MATRIX ELEMENT METHOD 

In this section, the measurement of the top quark mass 
using the Matrix Element method is described. The 
method is similar to the one of [5]; however, the calcula- 
tion of the signal probability has been revised, the nor- 
malization of the background probability is determined 
differently, and the method now allows a simultaneous 
measurement of the top quark mass and the jet energy 
scale [6]. 

An overview of the Matrix Element method is given 
in Section IV A. Section IV B describes how 6-tagging 
information is used in the analysis and discusses the pa- 
rameterization of the detector response. Details on the 
computation of the probabilities Psig and Pbkg are given 
in Sections IV C and IV D. 



A. The Event Probability 

To make maximal use of the kinematic information on 
the top quark mass contained in the event sample (or 
each individual event category in the case of the analysis 
using 6-tagging information), for each selected event a 
probability Pcvt that this event is observed is calculated 
as a function of the assumed top quark mass and jet en- 
ergy scale. The probabilities from all events are then 
combined to obtain the sample probability as a function 
of assumed mass and jet energy scale, and the top quark 
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mass measurement is extracted from this sample prob- 
ability. To make the probability calculation tractable, 
simplifying assumptions in the description of the physics 
processes and the detector response are introduced as 
described in this section. Before applying it to the data, 
the measurement technique is calibrated using fully sim- 
ulated events in the DO detector, and the assumptions in 
the description of the physics processes are accounted for 
by systematic uncertainties. 

It is assumed that the physics processes that can lead 
to the observed event do not interfere. The probability 
Povt then in principle has to be composed from probabil- 
ities for all these processes as 



(4) 



processes z 



where Pi is the probability for a given process i and fi 
denotes the fraction of events from that process in the 
event sample. In this analysis, Pevt is composed from 
probabilities for two processes, tt production and T4^-|-jets 
events, as 

Pevt{x; mtop,JES, /top) = /top Psig{x; mtop, JES) 

+ (l-/top)Pbkg(x; JP5) . (5) 

Here, x denotes the kinematic variables of the event, /top 
is the signal fraction of the event sample, and Pgig and 
Pbkg are the probabilities for ti and VK+jets production, 
respectively. The largest background contribution is from 
W-hjets events. Therefore, Pbkg is taken to be the prob- 
ability for VF+jets production. Contributions from QCD 
multijet events are not treated explicitly and are consid- 
ered as a systematic uncertainty. The signal probability 
Psig accounts for both possible flavor compositions in the 

hadronic W decay in £-|-jets tt events, W ud' and 
W cs': 



pW^ud I pW- 

2 v^sig ^^'^ 



(6) 



Because the event kinematics are the same, both final 
states are treated simultaneously in the probability cal- 
culation. 

To evaluate the tt probability, all configurations of tt 
decay products that could have led to the observed event 
X are considered. This includes different hadronic W de- 
cays as discussed above and all possible configurations y 
of the final state particles' four-momenta. The probabil- 
ity density for given partonic final state four-momenta 
y to be produced in the hard scattering process is pro- 
portional to the differential cross section da of the corre- 
sponding process, given by 



da{qq^tt^y;mtop) = 



{2nf\J^{qq~^tt^y)\^ 



d$fi 



(7) 



The symbol ^ denotes the matrix element for the pro- 
cess qq ^ ti ^ b{£i^)b{qq'), s is the center-of-mass en- 
ergy squared, qi and q2 are the momentum fractions of 



the colliding partons (which are assumed to be massless) 
within the colliding proton and antiproton, and dc&g is 
an element of six-body phase space. Here the symbol qq' 
stands for ud or cs'. 

To obtain the differential cross section d(j{pp tt ^ 
y, TTitop) in pp collisions, the differential cross section 
from equation (7) is convoluted with the parton density 
functions (PDF) for all possible fiavor compositions of 
the colliding quark and antiquark. 



da{pp ^tt^y; Wtop) (8) 
= y" E dqidq2 f{qi) f{q2) 



11, 12 

flavors 



da{qq tt ^ y; mtop) 



where f{q) denotes the probability density to find a par- 
ton of given flavor and momentum fraction q in the pro- 
ton or antiproton. 

The finite detector resolution is taken into account via 
a convolution with a transfer function W{x, y; JES) that 
describes the probability to reconstruct a partonic final 
state y as a; in the detector. The differential cross section 
to observe a given reconstructed ti event then becomes 



da{pp tt ^ x; mtop, JES) (9) 
= j da{pp -^ti^y; mtop)W{x, y; JES) . 



Because Pgig describes ^-|-jets tt events with both W 

ud' and W cs' decays, the transfer function depends 
on the quark flavors produced in the hadronic W decay 
when ^-tagging information is used. This is further dis- 
cussed in Section IV B. 

Only events that are inside the detector acceptance 
and pass the trigger conditions and offline event selec- 
tion are used in the measurement. The corresponding 
overall detector efficiency depends both on mtop and on 
the jet energy scale. This is taken into account in the 
cross section cjobs of ti events observed in the detector: 



o'obs(pp tt; mtop, JES) (10) 

= j da{pp ^ti^y; mtop)W{x,y; JES)f^c{x) , 
x,y 

where /acc ~ 1 for selected events and /acc = otherwise. 
The differential probability to observe a ti event as x 
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in the detector is then given by 

Psigix; mtop, JES) 

d(j{pp ^ tt ^ x; rritop, JES) 
cTobsipp tU m,op, JES) 
_ 1 

<^ohs{pP tt; mtop, JES) 

X / ^ dgidg2/(?i)/(92) 

91 > 92 



(11) 



flavors 



{2Tr)^ \^{qq^U^y)f 

W{x, y: JES) . 



d$6 



The parametrization of the matrix element and the com- 
putation of Pgig are described in Section IV C. 

Similarly, the differential background probability is 

computed as 



Pbkg(x; JES) 



(12) 



1 



o-obs(pp^ W+jets; JES) 



/ S 



dgidg2/(gi)/(g'2) 



9l>92 
flavors 



(27r)^ ^ W+jets 



91925 



■d$6 



where the matrix element and the total observed cross 
section for the process pp ^ W + jets have been used 
accordingly. Since the matrix element for VK+jets pro- 
duction does not depend on mtop, flakg is independent of 
mtop; however, Pbkg in principle does depend on the jet 
energy scale through the transfer function. Details about 
the Pbkg calculation can be found in Section IV D. 

To extract the top quark mass from a set of N mea- 
sured events xi,..,xn, a likelihood function is built from 
the individual event probabilities calculated according to 
Equation (5) as 



L{xi,..,xn] rutop, JES, f top) 



(13) 



JV 



= PevtiXi] mtop, JES, /top) . 



i=l 



For every assumed pair of values (mtop, JES), the value 
ftop^ that maximizes the likelihood is determined. The 
top quark mass and jet energy scale are then obtained by 
maximizing the likelihood 

L [xi, .., xn; mtop, JES, A'^ofKop, JES)) (14) 
= n ^-t {x^; mtop, JES, f^^^Ti^^^op, JES)) 



i=l 



with respect to mtop and JES, taking the correlation 
between both parameters into account. 



B. Description of the Detector Response 

The transfer function W{x, y; JES) relates the charac- 
teristics y of the final state partons to the measurements 
X in the detector. The symbol x denotes measurements 
of the jet and charged lepton energies or momenta and 
directions as well as ^-tagging information for the jets. 
A parameterization of the detector resolution is used in 
the probability calculation because the full GEANT-based 
simulation would be too slow. The full simulation is how- 
ever used to generate the simulated events with which the 
method is calibrated. 

The transfer function is assumed to factorize into con- 
tributions from each measured final state particle. The 
angles of all measured tt decay products as well as the 
energy of electrons are assumed to be well-measured; in 
other words, the transfer functions for these quantities 
are given by (5-distributions. This allows reducing the 
dimensionality of the integration over 6-particle phase 
space as described in Sections IV C and IV D. Conse- 
quently, contributions to the integral only arise if the 
directions of the quark momenta in the final state agree 
with the measured jet directions. In addition to the en- 
ergy resolution, one has to take into account the fact 
that the jets in the detector cannot be assigned unam- 
biguously to a specific parton from the tt decay. Conse- 
quently, all 24 permutations of jet-quark assignments are 
considered. 

In this section, the general form of the transfer func- 
tion in the topological and 6-tagging analyses is first dis- 
cussed, followed by a description of the jet energy and 
muon transverse momentum resolutions. 



1. Transfer Function in the Topological Analysis 

If no 6-tagging information is used, the transfer func- 
tion W{x, y; JES) is given by 



(15) 



W{x,y; JES) 
= W,({q/pT)TMIPT)T) 

24 4 
X — ^(5(angles) M^jet(i^jetj, ^'quarkfe; JES) , 



where and VFjot stand for factors describing the muon 
transverse momentum and jet energy resolutions, respec- 
tively. The sum is over the 24 different assignments of 
jets j to partons k. The factor (5(angles) denotes the 5 
distributions that ensure that assumed and reconstructed 
particle directions are identical, as discussed above. For 
e-|-jets events, the factor for the muon transverse momen- 
tum resolution is replaced with another ^-distribution. 
The neutrino is not measured in the detector and does 
not enter the transfer function. The jet transfer functions 
for light quark and charm jets are taken to be identical, 
and in the calculation of the background probability, all 
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jets arc assumed to be described by the light quark trans- 
fer function. Because the matrix elements for W ud 
and W — > cs' decays are equal, and the different flavor 
contributions to the ly+jets process are all parameter- 
ized by the VF-l-jets matrix element without heavy flavor 
quarks in the flnal state, no distinction between different 
processes is necessary in the topological analysis. 



2. Transfer Function in the b-tagging Analysis 

In the topological analysis, the information from the 

reconstructed jet energies determines the relative weight 
of different jet-parton assignments for a given partonic 
final state. The inclusion of 6-tagging information allows 
an improved identification of the correct permutation. 
This additional information enters the probability calcu- 
lation by weighting different permutations i of jet-parton 
assignments with weights Wi according to which jets, if 
any, are 6-tagged. This allows to give those permutations 
a larger weight that assign tagged jets to b quarks and 
untagged ones to light quarks. The transfer function is 
thus 
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FIG. 2: Monte Carlo study of the effect of charm-jet tagging 
on the signal to background probability ratio in the 6-tagging 
analysis, for tt events generated with mtop = 175 GeV that 
contain two ^-tagged jets. The Psig values arc calculated for 
the assumption mtop ~ 175 GeV. (a) Only the two jet-parton 
assignments in which tagged jets are assigned to b quarks are 
considered, (b) All weighted jet parton-assignments enter the 
probability calculation. In both plots, the hatched histogram 
corresponds to those cases where the two 6-tagged jets are 
correctly assigned to b quarks, which happens 84% of the 
time in the double tag sample. 



W{x,y;JES) (16) 

= w,({q/PTy;\{q/PT)r) 

24 4 

^U;i (5(angles) JJ Wjet(i0etj, Equarkk; JES) 
^ 24 

i=l 

(with a (5-distribution instead of the factor in the case 
of e+jets events). 

The weight Wi for a permutation i is parameterized as 
a product of individual weights Wj(j) for each jet. The 
latter are a function of the jet flavor hypothesis ak and 
the jet transverse energy Et.j and pseudorapidity rjj. For 
tagged jets, Wi{j) is equal to the per-jet tagging efficiency 
ejet(afc; Etj, T]j) where ak labels the three possible par- 
ton assignments to the jet: (a) h quarks, (b) c quarks, 
and (c) light quarks or gluons. For untagged jets, the 
Wi{j) factors are equal to 1 — ejet(Q;fc; Etj, rjj). If an 
event does not contain any 6-tagged jet, all the weights 
Wi{j) are set to 1.0. 

To compute the signal probability of events containing 
6-tagged jets, assumptions on the jet flavors are made for 
the calculation of the Wi such that hadronic W decays to 
ud and cs' final states need not be distinguished in the 
matrix element, allowing for a reduction of the computa- 
tion time. If an event contains exactly one 6-tagged jet, 
the quarks from the hadronic W decay are both assumed 
to be light quarks [u, d, or s). This is justified since the 
tagging efficiencies for 6 jets are much larger than those 
for other fiavors, and there arc two 6 jets per event. For 
events with two or more 6-tagged jets, a charm jet from 



the hadronic W decay is tagged in a non-negligible frac- 
tion of cases. Consequently, the quarks from the hadronic 
W decay are assumed to be charm quarks if the corre- 
sponding jet has been tagged, and light quarks otherwise. 

The need to include jet-parton assignments with a 
tagged charm jet in the probability calculation can be 
seen by comparing the signal and background probabili- 
ties. Figure 2(a) shows the ratio of signal to background 
probabilities calculated in a large sample of simulated 
ti events with two 6-tagged jets when only the two jet- 
parton assignments in which tagged jets are assigned to 
6 quarks are considered in the signal probability calcu- 
lation. The hatched histogram shows the correct assign- 
ments only, whereas the open histogram shows all combi- 
nations, including the ones in which a charm quark from 
the W decay was tagged. Figure 2(b) shows the same 
ratio when all combinations are included with their cor- 
responding weight as discussed above. The tail for low 
signal to background probability ratios in Fig. 2(a) arises 
because the correct jet-parton assignment is not included 
in the calculation in events where one of the tagged jets 
comes from a charm quark. It clearly shows the need to 
include these assignments in the signal probability calcu- 
lation. 

Because the different flavor contributions to the 
VF-l-jets process are parameterized by the W^-l-jets matrix 
element without heavy flavor quarks in the final state, the 
weights Wi for the background probability are all equal 
even if 6-tagged jets are present. Therefore, the back- 
ground probability calculated for the topological analysis 
is used in the 6-tagging analysis without modifications. 
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3. Parameterization of the Jet Energy Resolution 



transfer function for calorimeter jets, 
Eq ; JES) , yields the probability for a measure- 



The 

W'jot(^j 

ment Ej in the detector if the true quark energy is Eq 
For the case JES = 1, it is parameterized as 



W;^^{Ej,Eq-JES=l) 



2-K{p2+PzPb) 



(17) 



m^E^^ ^^^f_ [iE,-Eq)-p4- 



The parameters pi are themselves functions of the quark 
energy, and are parameterized as Unear functions of the 
quark energy so that 



Pi = Ui+Eq- bi 



(18) 



with 03 set to 0. 

The parameters and bi are determined from simu- 
lated events, after all jet energy corrections have been 
applied. The parton and jet energies are fed to an un- 
binned likelihood fit that minimizes the of the fit to 
Equation (17) with respect to and bi. A different set of 
parameters is derived for each of four 77 regions: I77I < 0.5, 
0.5 < |?7| < 1.0, 1.0 < |?7| < 1.5, and 1.5 < |?7| < 2.5, and 
for three different quark varieties: light quarks {u, d, s, 
c), b quarks with a soft muon tag in the associated jet, 
and all other b quarks. A total of 120 parameters de- 
scribe the transfer function for all jets, and are given in 
Tables IV and V. The transfer function for light quarks 
in the region jryj < 0.5 is shown in Fig. 3. 

For JES ^ 1, the jet transfer function is modified as 
follows: 



Wj,t{Ej,Eq-JES) 



Wi,,{^,Eq;l) 
JES 



(19) 
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TABLE IV: Light quark transfer function parameters (ai in 
GeV). 
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TABLE V: 6 quark transfer function parameters for jets with- 
out a muon (top) and for jets containing a muon (bottom) (a, 
in GeV). 



4.. Parameterization of the Muon Momentum Resolution 

To describe the resolution of the central tracking cham- 
ber, the resolution of the charge divided by the transverse 
momentum of a particle is considered as a function of 
pscudorapidity. The muon transfer function is parame- 
terized as 



w,({qiPTr;\{qiPT)r) 

1 ( {q/ptYi^'' - {q/pt) 



(20) 



27r(7 



exp 




where q denotes the charge and px the transverse mo- 
mentum of a generated (gen) muon or its reconstructed 
(rec) track. The resolution 



for \ri\ < rio 
\/(^o + [c{\v\-Vo)f for \v\ > Vo 



(21) 



is obtained from muon tracks in simulated events with 
the following values: 



(70 = 2.760 X 10"VGeV 

c = 5.93 X lO-^GeV 
r]o = 1.277. 



(22) 
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FIG. 3: Jet transfer functions for light quark jets, 0.0 < 
\ri\ < 0.5, for parton energies Ep = 30 GeV (solid), 60 GeV 
(dashed), and 90GcV (dash-dotted curve). The paramotriza- 
tion corresponds to the reference jet energy scale, JES = 1.0. 



(for the other case, replace b ^ b, £ ^ d, and v ^ u). 
Here, denotes the weak charge {Gp/V^ = gw/^^w)' 
rut and mw are the masses of the top quark (which is 
to be measured) and the W boson, and Tt and Tw are 
their widths. Invariant top and W masses in a particular 
event are denoted by nixyz and ruyz, respectively, where 
X, y, and z are the decay products. The cosine of the 
angle between particles x and y in the W rest frame is 
denoted by Cxy Hero and in the following, the symbols d 
and u stand for all possible decay products in a hadronic 
W decay. The top quark width is calculated as a function 
of the top quark mass according to [15]. 

The correct association of reconstructed jets with the 
final state quarks in Equations (24) and (25) is not 
known. Therefore, the transfer function takes into ac- 
count all 24 jet-parton assignments as described in Sec- 
tion IV B. However, in the case of the signal probabil- 
ity, the mean value of the two assignments with the 4- 
momenta of the quarks from the hadronic W decay inter- 
changed is computed explicitly by using the symmetrized 
formula 



The muon charge is not used in the calculation of fsig 
and Pbkg! however, for muons with large transverse mo- 
mentum it is important to take the possibility of charge 
misidentification into account in the transfer function. 



C. Calculation of the Signal Probability Psig 

The leading order matrix element for the process qq 
tt is taken to compute Psig- Neglecting spin correlations, 
the matrix element is given by [14] 



(23) 



where g^/{4:n) = is the strong coupling constant, /3 is 
the velocity of the top quarks in the tt rest frame, and Sgt 
denotes the sine of the angle between the incoming par- 
ton and the outgoing top quark in the tt rest frame. If the 
top quark decay products include the leptonically decay- 
ing W, while the antitop decay includes the hadronically 
decaying W, one has 



F = 



9u 



(1 



be 



) + mj^ (1 + cuf 



F = 



{mj^ - m^) -I- {niw^w) 



(24) 



(25) 




rn^f + {mw'i^wY 



(26) 



instead of (25), where only the terms containing are 
affected. Consequently, only a summation over 12 differ- 
ent jet-quark assignments remains to be evaluated. 

The computation of the signal probability Pgig involves 
an integral over the momenta of the colliding partons and 
over 6-body phase space to cover all possible partonic 
final states, cf. Equation (11). The number of dimensions 
of the integration is reduced by the following conditions; 

• The transverse momentum of the colliding par- 
tons is assumed to be zero. Conservation of 4- 

momcntum then implies zero transverse momen- 
tum of the tt system because the leading order ma- 
trix element is used to describe ti production. Also, 

the z momcntiim and energy of the tl system are 
known from the momenta of the colliding partons. 

• The directions of the quarks and the charged lep- 
ton in the final state are assumed to be exactly 
measured. 

• The energy of electrons from W decay is assumed 
to be perfectly measured. The corresponding state- 
ment is not necessarily true for high momentum 
muons, and an integration over the muon momen- 
tum is performed. 

After these considerations, an integration over the 
quark momenta, the charged lepton momentum (/x-|-jets 
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events only) , and the longitudinal component of the neu- 
trino momentum remains to be calculated. This calcu- 
lation is performed numerically with the Monte Carlo 
program VEGAS [16, 17]. The algorithm works most ef- 
ficiently if the one-dimensional projections of the inte- 
grand onto the individual integration variables have well- 
localized peaks. The Breit-Wigner peaks of the integrand 
corresponding to the two top quark and two W boson de- 
cays in the tt matrix element are more localized than the 
peaks from the jet transfer functions, suggesting that the 
masses are better integration variables leading to faster 
convergence. The computation of the parton kinemat- 
ics from the integration variables must however be per- 
formed in each integration step, and this task simplifies 
to solving a quadratic equation when choosing Pz,bi' as 
an integration variable instead of the mass of the lep- 
tonically decaying W (both solutions of the quadratic 
equation are considered when determining Pgig)- There- 
fore, the following integration variables are chosen for the 
computation of Pgig: 

• the magnitude \pd \ of the momentum of one of the 
quarks from the hadronic W decay, with < \pd\ < 
500 GeV, 

• the squared mass of the hadronically decaying 
W,0<m%.< (400 GeV)2, 

• the squared mass nT^^- of the top quark with the 
hadronic W decay, < < (500 GeV)^, 

• the squared mass ml^^ of the top quark with the 
leptonic W decay, < mg^^ < (500 GeV)^, 

• the z component Pz,bv of the sum of the momenta 
of the b quark and neutrino from the top quark 
with the leptonic W decay, —500 GeV < Pz,bv < 
-h500 GeV, and 

• the muon charge divided by the muon trans- 
verse momentum (in the /z-|-jets channel only), 
-1/(100 MeV) < (g/pr)^ < +1/(100 MeV). 

Thus, for each point in the (|pd|, 171%^, m?^_, mg^^, 
Pz.bi^ [, ('z/pt)^ ]) integration space the following compu- 
tation is performed for each of the 12 possible jet-parton 
assignments (where the symmetrized form of the matrix 
element according to Equation (26) is used): 

1. The 4- momenta of the tt decay products arc calcu- 
lated from the values of the integration variables, 
the measured jet and lepton angles, and the elec- 
tron energy (in the e+jets case). 

2. The matrix element is evaluated according to Equa- 
tions (23), (24), and (26). 

3. The parton distribution functions are evaluated. 
For consistency with the leading-order matrix ele- 
ment, we use the CTEQ5L [18] parton distribution 
functions, summing over all possible quark flavors. 



4. The probabilities to observe the measured jet en- 
ergies and muon transverse momentum given the 
energies and momentum computed in the first step 
are evaluated using transfer functions. 

5. The Jacobian determinant for the transformation 
from momenta in Cartesian coordinates to the {\pd\, 
m^, m2^_, ml^^, Pz^ [> (^/Pt)^ ]) integration 
space is included. 

The precision of the Pgig calculation varies from typically 
2 % to a maximum of 10 %. 

To normalize the signal probability, the integral 
/ (i(j{pp tt ^ x; rritop, J ES)f s,cc{x) over 16-dimen- 
sional phase space has been computed as a function of 
mtop and JES. The detector acceptance and efficiency 
is taken into account as outlined in Equation (10). The 
results are shown in Fig. 4 for e+jets and /i+jets events 
as a function of mtop for various choices of the JES scale 
factor. 



D. Calculation of the Background Probability Pbkg 

To calculate Phkg, the jet directions and the charged 
electron or muon arc taken as well-measured. The in- 
tegral over the quark energies in Equation (12) is per- 
formed by generating Monte Carlo events with parton 
energies distributed according to the jet transfer func- 
tion. In these Monte Carlo events, the neutrino trans- 
verse momentum is given by the condition that the trans- 
verse momentum of the M^+jets system be zero, while 
the invariant mass of the charged lepton and neutrino is 
assumed to be equal to the W mass to obtain the neu- 
trino z momentum (both solutions are considered). The 
VECBOS [19] parameterization of the matrix clement is 
used. The mean result from all 24 possible assignments 
of jets to quarks in the matrix element is calculated. A 
minimum of 10 Monte Carlo events is generated for each 
measured event x, and the relative spread of the resulting 
Pbkg values is evaluated as the standard deviation divided 
by the mean. If the relative spread is larger than 10%, 
another 10 Monte Carlo events are evaluated, and this 
procedure is repeated until a 10% relative uncertainty is 
reached or a maximum of 100 Monte Carlo events has 
been considered. 

To normalize the background probability density, 
<^ohs{pp ^ W + jets; JES) is chosen such that the to- 
tal signal fraction /top in the analysis without b tagging 
is reproduced in the fit to simulated event samples con- 
taining ti and W+jets events. This makes use of the fact 
that /top is underestimated in the fit if the background 
probabilities are too large and vice versa. 

In the simulation, about 20 — 30 % of tt events have 
jets and partons that cannot be unambiguously matched, 
i.e., at least one of the four reconstructed jets cannot be 
matched to a parton from the tt decay within ATZ < 0.5. 
These events yield poor top mass information and de- 
grade the uncertainty estimate of the likelihood fit. Fig- 
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FIG. 4: Observed tt cross section computed with the leading- 
order matrix element for (a) e+jets and (b) n+jets events as 
a function of the top quark mass mtop for different choices of 
the JES scale factor: JES = 1.12 (dash-dotted), JES = 1.0 
(solid), and JES = 0.88 (dotted lines). 



urc 5 illustrates that jct-parton matched events tend to 
have a higher signal than background probability density, 
which is how the mass fit identifies them as signal-like. 
There is no such separation for signal events in which 
one or more jets cannot be matched to a parton, so that 
these events contribute much less mass information to 
the final likelihood. This observation is consistent with 
the fact that a leading-order matrix element is used to 
describe tt events. Therefore, only jet-parton matched 
events are used to calibrate the -Pbkg normalization. On 
average, the /top fit will consequently yield the fraction of 
jet-parton matched (leading-order) tt events in the event 
sample. The quoted /top values are corrected for this 
effect. 

The Pbkg normalization is determined as follows: 

• A large ensemble of simulated tt and VF+jets events 
is composed with the signal fraction as determined 



FIG. 5: Distributions of log]^Q(Psig/-Pbkg) for tt events with 
mtop = 175 GeV (solid) and W^-l-jets events (dashed lines) 
for (a) e-l-jets events and (b) /i-|-jets events. The Psig val- 
ues arc calculated for the assumption mtop ~ 175 GcV. The 
distributions for signal and background events are normalized 
individually. The distributions for those tt events that fail the 
requirement of jets matched to partons are shown separately 
(dash-dotted lines). 



by the topological likelihood fit described in Sec- 
tion IIIG. 

• The top mass likelihood fit described in Section V 

is applied to the sample and the -Pbkg normalization 
is adjusted iteratively until the fit result yields the 
true signal fraction. 

• The normalization of Pbkg cannot depend on the 

top quark mass. Therefore, the above steps are 
applied to ti Monte Carlo samples with different 
generated masses. The mean of all results is taken 
as the Pbkg normalization. 

This procedure is applied separately for e+jets and 
li+iets events. Note that the topological likelihood dis- 
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criminant is only used to determine the normalization 
of the background probability and the sample composi- 
tion for ensemble tests used to calibrate the procedure. 
The topological likelihood discriminant does not other- 
wise enter the top quark mass fit. 



V. TOP QUARK MASS MEASUREMENT 
USING TOPOLOGICAL INFORMATION 

A. Top Quark Mass Fit 

The top quark mass and overall jet energy scale JES 
are determined as optimal values of the likelihood for the 
sample of selected events, which depends on the Psig and 
Pbkg values. For each measured event, Pgig is calculated 
for various values of mtop in steps of 2.5 GeV and various 
values of JES in steps of 0.01. It has been found that 
it is not necessary to compute the background probabil- 
ity for different values of the jet energy scale. There- 
fore, all Pbkg values arc computed for JES = 1 only. 
Both Pgig and Pbkg are normalized as described in Sec- 
tions IV C and IV D, using separate constants for e-|-jets 
and y^-|-jets events. The top quark mass measurements on 
the e-|-jets, /K-|-jets, and combined ^-|-jets event samples 
are in each case derived from the likelihood of the event 
sample, given by Equation (13), in the way described 
below. 

For given values of rntop and JES, each event prob- 
ability Pevt = /top Psig + (1 - /top) Pbkg depends on the 
signal fraction /top of the sample, and consequently, the 
value of the likelihood for the event sample is a function 
of /top- For each (rntop,JES) parameter pair, the best 
/top parameter value is determined, and the likelihood 
value corresponding to this value is used in further com- 
putations. The overall result quoted for the fitted signal 
fraction /top is derived from the value obtained at the 
point in the grid of {mtop,JES) assumptions with the 
maximum likelihood value for the event sample. The un- 
certainty on /top is compiited by varying /top at fixed 
mtop and JES until A(— InL) = This uncertainty 
does not account for correlations between /top, mtop, and 
JES. 

The result for the top quark mass is obtained from a 
projection of the two-dimensional grid of likelihood val- 
ues onto the mtop axis. In this projection, the correlation 
between mtop and the JES parameter is taken into ac- 
count. The probability for a given mtop hypothesis is ob- 
tained as the integral over the likelihood as a function of 
JES, using linear interpolation between the grid points 
and Gaussian extrapolation to account for the tails for 
JES values outside the range considered in the grid. 

The probabilities as a function of assumed top mass are 
converted to — InL values. These — InL points are then 
fitted with a fourth order polynomial in the region defined 
by the condition A In i < 3 around the best value. The 
— In L points on either side of the A In i < 3 region arc 
each fitted with a parabola, and Gaussian extrapolation 



is used to describe the tails outside the range of mtop 
hypotheses considered. The mtop value that maximizes 
the fitted probability is taken to be the measured value of 
the top quark mass. The lower and upper uncertainties 
on the top mass are defined such that 68% of the total 
probability integral is enclosed by the corresponding top 
mass values, with equal probabilities at both limits of the 
68% confidence level region. 

The same projection and fitting procedure is applied 
to determine the value of the JES parameter. 



B. Validation of the Method 

The method is first validated using parton-level simu- 
lated tt and M^-hjets events. These have been generated 
with leading-order event generators (madgraph [20] for 
<< events, ALPGEN for VK+jets events), i.e., no initial or 
final state radiation is included. The jet energies in these 
events are smeared according to the transfer functions 
described in Section IV B (the treatment of the muon 
transverse momentum integration has been checked with 
additional ensemble tests not described here). 

Ensembles are composed with 75 events, 40 % of which 
are tt signal events. A total of 1000 events for top masses 
of 160, 170, 175, 180, and 190 GeV each are used, along 
with 1000 M^-hjets events. In addition, samples with 
TOtop = 175 GeV with all jet energies scaled by 0.95 and 
1.05 are prepared in order to validate the JES fit result. 
All events are required to pass the kinematic selection cri- 
teria listed in Table I. The final state jets and the charged 
lepton must be separated according to ATZ{j,j') > 0.5 
and A7?.(£,j) > 0.5. The signal normalization is obtained 
according to this selection, see Section IV C. mtop and 
JES are obtained for each ensemble as described in Sec- 
tion VA. The results of this test show that the fitted 
top mass and jet energy scale are unbiased within sta- 
tistical uncertainties of 300 MeV and 0.003, respectively. 
Furthermore, the fitted mtop value does not depend on 
the input JES value used in the ensemble generation, 
and similarly, the fitted JES value is independent of the 
true input top mass. 

To test that the uncertainties obtained from the fit de- 
scribe the actual measurement uncertainty, the deviation 
of the fitted top mass from the true value is divided by 
the fitted measurement uncertainty. The upper (lower) 
measurement uncertainty is taken if the fitted value is 
below (above) the true value. This definition is chosen 
to account for the possibility of asymmetric uncertain- 
ties. This distribution of deviations normalized by the 
measurement uncertainty is fitted with a Gaussian, and 
its width, commonly referred to as "pull width," is in 
agreement with 1.0. This is also the case for the jet en- 
ergy scale measurement, for which same test has been 
performed. 

The events used in the test outlined in this section 

have been generated with the same simplified model that 
is used in the probability calculation for the description 
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of the production of signal and background events and 
the detector response, cf. Section IV. As it cannot be 
assumed that this simphfied model correctly reproduces 
every aspect of the data, the method for measuring the 
top quark mass has been calibrated with Monte Carlo 
events that have been generated with the full DO simula- 
tion. Any deviations observed in this calibration step are 
taken into account in the final result. The calibration is 
described in the following section. 
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C. Calibration of the Method 

The default DO Monte Carlo events, generated as de- 
scribed in Section III and passed through the full simu- 
lation of the DO detector, are found to describe the data 
well. They are therefore used to derive the final cal- 
ibration of the fitting procedure, tt samples with top 
quark masses of 160, 170, 175, 180, and 190 GeV and 
a W^-f-jets sample are used. In addition, samples with 
'T^top = 175 GeV where all jet energies are scaled by 0.92, 
0.96, 1.04, and 1.08 are prepared in order to calibrate 
the JES fit. For each sample and each lepton chan- 
nel (e-|-jets and /x-|-jets) , Pgig and Pbkg are calculated for 
1000 events which pass the event selection. Ensembles 
arc drawn from these event pools, with an ensemble com- 
position as measured for the data sample. Each proba- 
bility is normalized according to the flavor of the isolated 
lepton (see Sections IV C and IV D). The QCD contri- 
bution is not added during the calibration but treated as 
a systematic uncertainty (cf. Section VII). 

In Fig. 6, 68% confidence interval distributions are 
shown for ensembles with rutop = 175 GeV and JES = 
1.0. For many pseudo-experiments, the true mtop (JES) 
value is expected to be within the fitted uncertainties 
in 68% of the pseudo-experiments (corresponding to a 
value of the conficcncc interval distribution of 0.68 for 
this value), while other mtop (JES) values should be less 
likely to be within the uncertainties. When the error 
interval resulting from the integration of 68% of the like- 
lihood distribution does not include the input top mass 
{JES) value 68% of the time, the uncertainty is inflated 
to correspond to an integration over a larger interval. 
The calibration results for the combined fit to the e-|-jets 
and /x-|-jets ensembles arc shown in Figs. 7 and 8. The 
fit results are corrected for the offsets o and slopes s, and 
for the deviations of the pull width w from 1.0 given in 
Table VI to obtain the final results and their statistical 
uncertainties as follows: 



'Tltop 

Amtop 
JES 
AJES 



mf^p - On 



- 175 GeV 



+ 175 GeV , 



(Amtop) 

JES^' - ojES - 1 



SJES 



+ 1 , and 



wjES {AJES) 



fit 



(27) 



FIG. 6: Calibration of the Matrix Element mass fitting proce- 
dure for the topological analysis, using ensembles with a top 
quark mass of 175 GeV and JES = 1.0. The 68% confidence 
interval distributions for (a) the measured top quark mass 
and (b) the jet energy scale is given by the solid, the upper 
and lower error bands by the dashed histograms. A value of 
0.68 as indicated by the dash-dotted line would mean that 
the corresponding mtop {JES) value is included in the fitted 
68% confidence interval in 68% of the ensembles. 





offset 


slope s 


pull width w 


mtop 


1.375 ± 0.085 GeV 


1.034 ± 0.011 


1.06 ± 0.01 


JES 


-0.028 ± 0.001 


0.934 ± 0.021 


1.09 ± 0.01 



TABLE VI: Calibration of the Matrix Element mass fitting 
procedure for the topological analysis. The offsets are quoted 
for a true top quark mass of 175 GeV and a true jet energy 
scale of 1.0, respectively. Only statistical uncertainties are 
quoted in this table. 



D. Result 



The Matrix Element method is applied to the 370 pb~^ 
lopton+jets data set. The calibrations for m,top derived 
in the previous section are taken into account. The cal- 
ibrated fit result for the combined lepton-l-jets sample 
is shown in Fig. 9. In this figure, the probability as a 
function of assumed top mass is shown together with the 
fitted curve (the polynomial fitted to the — In L values as 
described in Section VA has been transformed accord- 
ingly), and the central value and 68% confidence level 
interval are indicated. The probability as a function of 
assumed JES parameter is also shown. The top quark 
mass is measured to be 



e-l-jets 
'''top 

'''■top 

^-l-jets 



m. 



top 



167.0 
173.0 
169.2 



(stat. + JES) GeV 



+ 7.3 
-10.6 



"^")-J (stat. + JES) GeV 



-10. i 

+5.0 
-7.4 



(stat. + JES) GeV . (28) 



The statistical uncertainties are consistent with the ex- 
pectation. A comparison of the fitted uncertainties on 
mtop and JES with the expectations from ensemble tests 
is given in Fig. 10. The fit yields a signal fraction /top of 
0.33 ± 0.06 (stat.), in good agreement with the result of 
the topological likelihood fit. The fitted jet energy scale 
is JES = 1.0481o;q4q (stat.) and indicates that the data 
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FIG. 7: Calibration of the Matrix Element mass fitting pro- 
cedure for the topological analysis. The upper plots show the 
reconstructed top mass (a) and the measured jet energy scale 
(b) as a function of the input top mass. The two lower plots 
show the reconstructed top mass (c) and the measured jet 
energy scale (d) as a function of the input jet energy scale. 
The solid lines show the results of linear fits to the points, 
which are used to calibrate tlie measurement technique. The 
dashed lines would be obtained for equal fitted and true values 
of mtop and JES. 
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FIG. 8: Calibration of the Matrix Element mass fitting pro- 
cedure for the topological analysis. The upper plots show the 
widths of the pull distributions for the top mass (a) and jet 
energy scale (b) as a function of the input top mass. The two 
lower plots show the widths of the pull distributions for the 
top mass (c) and jet energy scale (d) as a function of the input 
jet energy scale. The solid lines show the mean pull width, 
while the dashed lines indicate a pull width of 1.0. 



is consistent with the reference scale. 

For a fixed jet energy scale, the statistical uncertainty 
of the fit is 2 GeV; thus the component from the jet 
energy scale uncertainty is lg'7 GeV. Systematic uncer- 
tainties are discussed in Section VII. 

To show the likelihood as a function of both rutop and 
JES simultaneously, the — In L values have been fitted 
with a two-dimensional fourth-degree polynomial with its 
minimum fixed to the measurements mentioned above. 
The resulting contours corresponding to Aln_L = 0.5, 
2.0, 4.5, and 8.0 relative to the minimiun are shown in 
Fig. 11. Note that the statistical measurement uncer- 
tainties quoted on mtop and JES are obtained from the 
one-dimensional projections as discussed above; Fig. 11 
therefore serves only illustrative purposes. Because of 
non-Gaussian tails, the projections of the AlnL = 0.5 
contour shown in Fig. 11 onto the mtop and JES axes 
do not exactly correspond to these quoted statistical un- 
certainties. 



VI. TOP QUARK MASS MEASUREMENT 
USING B JET IDENTIFICATION 



A. Top Quark Mass Fit 

The incorporation of 6-tagging information introduces 
two significant modifications to the Matrix Element mass 
fitting technique. First, 6-tagging information is used 
to determine the relative weights Wi of the difii'erent jet- 
parton assignments in the signal probability calculation. 
The Wi are parameterized as a function of the jet trans- 
verse energy Et and pseudorapidity r] and the assumed 
flavor ak of the parton corresponding to the jet, as de- 
scribed in Section IV B 2. The signal probability is then 
computed according to Equation (11). The background 
probability is identical to that used in the topological 
analysis according to Equation (12). 

The second modification is to classify events into three 
categories according to the number of 6-tagged jets. Each 
of these categories will have different signal fractions and 
background compositions due to the relative suppression 
of W-|-jets events with dominantly light quark and gluon 
jets. The event categories are exclusive and correspond 
to i) no &-taggcd jet, ii) exactly one tagged jet, and iii) 
two or more tagged jets. 

When the analysis is separately performed in each ntag 
category, the signal fractions /t"'p* are determined inde- 
pendently for each category, and P^vt'^ is calculated as 

PZX'^ix; mtop, JES, /top) = ^ P^^' i^; ^top, JES) 
+{l-f::;^)P^^^ix;JES). (29) 

To combine the three categories into one analysis the 

three purities /t"*p'' have to be related to one inclusive 
signal purity /top . The purity of the ntag sample is given 
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FIG. 9: Application of the topological Matrix Element 
method to the data. The mtop and JES axes correspond 
to the calibrated values. Plot (a) shows the probability as a 

function of assumed top quark mass. The correlation with 
the jet energy scale is taken into account. The fitted curve is 
shown, as well as the most likely value and the 68% confidence 
level region. The corresponding plot for the JES parameter 
is shown in (b). 



FIG. 10: Uncertainties on mtop (a) and JES (b) in the topo- 
logical analysis. The distributions of fitted uncertainties ob- 
tained from ensemble tests are shown by the histograms. Both 
upper and lower uncertainties are shown; their distributions 
are very similar. The upper (lower) uncertainty in the data is 
indicated by the solid (dashed) arrow. The probability for a 
lower uncertainty on mtop with a magnitude larger than that 
observed in the data is 2%. 



by 



as 



top 



iy-"tag 

sig 



iy"tag I iy"tag 

sig ~ bkg 



(30) 



The numbers of signal and background events after b tag- 
ging, A^^g** and N^^^ , can be related to the correspond- 
ing numbers for the inclusive sample, A^sig and A^bkg, by 



sig 



iVsi 



and iV^C " ^^'^K 



bkg 



(31) 



where e^;°« and e^l^^ 



arc the average tagging efficiencies 
for signal and background, respectively. They are defined 




''{{ai} = bbud) + e>M{ai} = bbcs)) 



(32) 



with relative fractions /$ of the different flavor contri- 
butions $ to the VF+jets background as given in Ta- 
ble III. The jets in selected QCD multijet background 
events have kinematic characteristics similar to those of 
jets in selected IF-l-jets background events. Concerning 
the event 6-tagging probabilities, we therefore do not 
distinguish between QCD multijet and W+jcts back- 
grounds. The difference between multijet and M^+jets 
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of event-tagging efficiencies and signal purity before tag- 
ging. In order to extract the top quark mass from the 
total sample of selected events, the likelihoods in the in- 
dividual event categories are then combined as 

L{xi, .., xn; mtop, JES, /top) (36) 

= n II p::r{xi;m,op, JES, , 

ntag=0,l,>2 i=l 

where A'^"**^ is the number of events in each of the three 

tag-categories. As in the topological analysis, we de- 
termine the value /to'jf* (n2top, JES) that maximizes the 
likelihood L in Equation (36) for each pair of assumed 
values of mtop and JES. The top quark mass and jet 
energy scale are then obtained by maximizing 



FIG. 11: Application of the topological Matrix Element 
method to the data. Fit of a two-dimensional fourth-order 

polynomial to the — In L values as a function of both mtop and 
JES. Shown are the contours corresponding to AlnL = 0.5, 
2.0, 4.5, and 8.0 relative to the minimum. 







Subsample 


Channel 


0-tag 


1-tag > 2-tag 


e-|-jets 


0.68 


0.86 0.94 


/^-|-jets 


0.77 


0.84 0.90 



TABLE VII: Signal purity correction factors c"*"* . 



kinematics is treated as a systematic uncertainty. The 
relation between /t"'p* and the inclusive signal purity 
/top = -^Vsig/(-^Vsig + iVbkg) is then 



top 



/top?- 



/top(r'^ 



1) + 1 



where 



"tag 
^"tag _ Sig 
' ^"tag 
^bkg 



(33) 



(34) 



Equation (33) needs to be corrected for the fact that the 
fraction of tt events that are jet-parton matched is differ- 
ent in each tag-multiplicity sample. Thus, a correction 
factor c"*'*i5 defined as the ratio of fitted tt fraction over 
expected tt fraction is introduced as an intercalibration 
of the /top" values. The top fraction for a given tag- 
multiplicity sample is then 



/■"tag /."tag 

Jtop — 



/ 



top' 



/top(r" 



1 



1 



(35) 



The correction factors c"""5 are different for e-|-jets and 
/x-|-jets events, and are given in Table VII. 

Equation (35) defines the dependence of the signal pu- 
rity on tagging multiplicity, as a function of the ratio 



L{xi,..,xn; TOtop, JES, /top '(mtop, JES)) 



(37) 



iV^tag 



= n IKvr m,op, JES, /t'j.f ""^ (mtop, JES)) 

"tag=0,l,>2 j=l 

as described in Section V A (for the 0-tag sample, the fit 
range is restricted to AlnL < 1). 



B. Calibration of the Method 

The calibration is obtained following a similar proce- 
dure as described in Section VC. The number of events 
in each tag-multiplicity class is calculated by multipliying 
the expected number of selected events (before tagging) 
by the average event-tagging probability in each process: 



^"tag 
A/-"tag 



N!r^;r ^^d 



ArSol "tag 

^^bkg^lV+jots ■ 



(38) 



The VF-fjets background is classified into two categories 
according to the differences in event kinematics: Wjjjj 
{W + 4 jets without heavy fiavor) and Whf {W + 4 jets 
including heavy flavor). Thus, the background composi- 
tion for the ensemble tests after tagging is given by 



AT"tag 
AT"tag 



= ^bkl /l 



bkg 
-'^bkg 



I Wjjjj ^wJjjj ^nd 



E . 



$ fid 



(39) 



where fwjjjj is the fraction of T4^ + light jets, $ denotes 
one of the five VF-l-heavy flavor subprocesses (see Ta- 
ble III), and /$ is the corresponding fraction of these 
subprocesses. Table VIII shows the ly+jcts composition 
in the 0, 1, and > 2 tag samples used in the ensemble 
tests. The average number of events in each tag-catc!gory 
and for each sample are fluctuated according to a Poisson 
distribution. For a given tag-category ritag, the decision 
of which jets are tagged is made by randomly selecting 
^tag jets as tagged jets, taking into account the Et and 
rj dependence of the tagging efficiencies. 
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Subsamplo 



Contribution 


0-tag 


1-tag 


> 2-tag 


Wjjjj 


90.9% 


19.4% 


0.0% 


Whf 


9.1% 


80.6% 


100.0% 



TABLE VIII: Background composition used in the ensem- 
ble tests for the 6-tagging analysis. The contribution from 
ly+jets events without heavy flavor is given in the first line, 
the contribution from events with heavy flavor jets in the sec- 
ond line. 





offset 


slope s 


pull width w 


mtop 


1.932 ± 0.085 GeV 


1.018 ±0.011 


1.11 ±0.01 


JES 


-0.028 ± 0.001 


0.945 ± 0.021 


1.09 ±0.01 



TABLE IX: Calibration of the Matrix Element mass fitting 
procedure for the ^-tagging analysis. The offsets are quoted 
for a true top quark mass of 175 GeV and a true jet energy 
scale of 1.0, respectively. Only statistical uncertainties are 
quoted in this table. 



In Fig. 12, 68% confidence interval distributions are 
shown for ensembles with mtop = 175 GeV and JES = 
1.0. The calibration results for the combined fit are 
shown in Figs. 13 and 14. The final fit results are cor- 
rected for the biases and for the deviation of the pull 
width from 1.0 given in Table IX. 
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FIG. 13: Calibration of the Matrix Element mass fitting pro- 
cedure for the 6-tagging analysis. The upper plots show the 
reconstructed top mass (a) and the measured jet energy scale 
(b) as a function of the input top mass. The two lower plots 
show the reconstructed top mass (c) and the measured jet 
energy scale (d) as a function of the input jet energy scale. 
The solid lines show the results of linear fits to the points, 
which are used to calibrate the mcasurcrncnt technique. The 
dashed lines would be obtained for equal fitted and true values 
of mtop and JES. 



C. Result 

The Matrix Element 6-tagging method is applied to 
the same event sample as in Section VD with a cali- 
bration according to the results from Section VI B. The 
probability is shown in Fig. 15 as a function of mtop and 
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FIG. 12: Calibration of the Matrix Element mass fitting pro- 
cedure for the 6-tagging analysis, using ensembles with a top 
quark mass of 175 GeV and JES = 1.0. The 68% confidence 
interval distributions for (a) the measured top quark mass 
and (b) the jet energy scale is given by the solid, the upper 
and lower error bands by the dashed histograms. A value of 
0.68 as indicated by the dash-dotted line would mean that 
the corresponding mtop {JES) value is included in the fitted 
68% confidence interval in 68% of the ensembles. 
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FIG. 14: Calibration of the Matrix Element mass fitting pro- 
cedure for the 6-tagging analysis. The upper plots show the 
widths of the pull distributions for the top mass (a) and jet 
energy scale (b) as a function of the input top mass. The two 
lower plots show the widths of the pull distributions for the 
top mass (c) and jet energy scale (d) as a function of the input 
jet energy scale. The solid lines show the mean pull width, 
while the dashed lines indicate a pull width of 1.0. 
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JES hypothesis for each of the three tag categories. The 
central values and the 68% confidence level intervals are 
indicated in the figures. The individual results for the 
top quark mass are 



top 

«+jets 
top 

mL+j^*^(2-tag) = 



TO, 



TO, 



(0-tag) = 
(1-tag) = 



174.4tl|^ 



173.1 
163.2 



+5.1 
-5.2 

+6.8 
-6.2 



"top 

in the 0-tag, 
corresponding 
JES{Q-i&g) = 
1.01lt°:°l^(stat 
respectively. 

The result for the combined event 
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(40) 



1-tag, and > 2-tag 
results for the jet 
= 0.986ll];°^t (stat.), 
, and J£'S'(2-tag) = 



categories. The 

energy scale are 
J£;S'( 1-tag) = 
l-094tn66(stat.), 



e+iets 
^top 

u+jots 
Hop 



TO. 



£+jets 
top 



170.lt^j(stat. 
172.6l^:e (stat. 
170.3t4;5 (stat. 



sample is 

- JES) GeV 
^ JES) GeV 

- JES) GeV; 



(41) 



the f+jets measurement is shown in Fig. 16. The sta- 
tistical uncertainties are consistent with the expectation. 
Fi gure 17 shows the distributions of the expected TO.top 
uncertainty compared to the observed result. The fit 
yields a signal fraction /top of 0.31 ± 0.09 (stat.), in good 
agreement with the result of the topological likelihood fit. 
The fitted jet energy scale is JES = 1.027lJJ;^^^ (stat.) 
and indicates that the data is consistent with the refer- 
ence scale. 

For a fixed jet energy scale, the statistical uncertainty 
of the fit is ±2.5 GeV; thus the component from the jet 
energy scale uncertainty is j^g'y GeV. Systematic uncer- 
tainties arc discussed in Section VII. 

To show the likelihood as a function of both mtop and 
JES simultaneously, the — In L values have been fitted 
with a two-dimensional fourth-degree polynomial with its 
minimum fixed to the measurements mentioned above. 
The resulting contours corresponding to AlnL = 0.5, 
2.0, 4.5, and 8.0 relative to the minimum arc shown in 
Fig. 18. Note that because of non-Gaussian tails, the pro- 
jections of the Alni = 0.5 contour onto the TOtop and 
JES axes do not exactly correspond to the 68% con- 
fidence intervals around the most likely values; Fig. 18 
therefore serves only illustrative purposes. 



VII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties arise from three sources: 
modeling of the physics processes for tt production and 
background, modeling of the detector performance, and 
uncertainties in the methods themselves. Table X lists 
all uncertainties. The jet energy scale uncertainty is in- 
cluded in the statistical uncertainty. The total systematic 
uncertainty on the top mass measurement is obtained 
by adding all contributions in quadrature. In general, 
to evaluate systematic uncertainties, the simulation of 
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FIG. 15: Application of the Matrix Element 6-tagging method 
to the data. The fitted mtop and JES likelihoods for each of 

the 3 tag-categories: 0-tag ((a) and (b)), 1-tag ((c) and (d)), 
and > 2-tag ((e) and (f)). The 68% confidence-level interval 
around the most likely value is shown by the hatched region 
under the fitted curve. 



events uscxl to calibrate the measurement has been var- 
ied, while the measurement method itself has been kept 
unchanged. 



A. Physics Modeling 

Signal modeling: When tt events arc produced 
in association with a jet, the additional jet can be 
misinterpreted as a product of the tt decay. Also, 
the tt system may then have significant transverse 
momentum, in contrast to the assumption made in 
the calculation of Pgig- In spite of the event selec- 
tion that requires exactly four jets, these events can 
be selected if one of the jets from the ti decay is not 
reconstructed. 

Such events arc present in the simulated events 
used for the calibration of the method. To assess 
the uncertainty in the modeling of these effects, 
events have been generated using a dedicated sim- 
ulation of the production of ti events together with 
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FIG. 16: Application of the Matrix Element 6-tagging method 
to the data. The final results of the fitted mtop (a) and JES 
(b) Ukelihoods for the combined event sample are shown. The 
68% confidcncc-lovcl interval around the most likely value is 
indicated by the hatched region under the fitted curve. 



FIG. 17: Uncertainties on mtop (a) and JES (b) obtained in 
the 6-tagging analysis with the combined sample. The distri- 
butions of fitted uncertainties obtained from ensemble tests 
arc shown by the histograms. Both upper and lower uncer- 
tainties arc shown; their distributions are very similar. The 
upper (lower) uncertainty in the data is indicated by the solid 
(dashed) arrow. 



an additional parton. The fraction of such events 
is estimated to be no larger than 30% (according 
to the difference between cross section calculations 
in leading and ncxt-to- leading order). 

Two large ensembles of simulated events are com- 
posed according to the sample composition in the 
data, one using only events with an additional par- 
ton for the signal, and the second with the default 
simulation. The result obtained with the default 
calibration is quoted as central value. A systematic 
uncertainty of 30% of the difference in top mass re- 
sults between these two ensembles is quoted. 

In addition, simulated gg tt and qq tt events 
have been compared. The top mass calibration has 
been rederived using only gg — > tt or qq — > events 
to simulate the signal, and no significant difference 
has been found. Thus no additional uncertainty on 



the result is assigned. 

• Background modeling: In order to study the 

sensitivity of the measurement to the choice of 
background model, the standard VF-l-jets Monte 
Carlo sample is replaced by an alternative sam- 
ple with the default factorization scale of = 
m'^Y + J^jP'r.] replaced by Q'^ = (ptj)'^- One large 
ensemble of events is composed using both the de- 
fault and the alternative background model. The 
difference of results obtained with these ensembles 
is symmetrized and is assigned as a systematic un- 
certainty. 

• PDF uncertainty: Leading-order matrix ele- 
ments are used to calculate both Pgig and Pbkg- 
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FIG. 18: Application of the Matrix Element ?>-tagging method 
to the data. Fit of a two-dimensional fourth-order polynomial 

to the — InL values as a function of both mtop and JES. 
Shown are the contours corresponding to AlnL = 0.5, 2.0, 
4.5, and 8.0 relative to the minimum. 



Source of Uncertainty 



Topological 
Analysis 



6-Tagging 
Analysis 



Statistical uncertainty 






+4.1 -4.5 


and jet energy scale 


-1-5.0 


-7.4 


Physics modeling: 








Signal modeling 


±0 


34 


±0.46 


Background modeling 


±0 


32 


±0.40 


PDF uncertainty 


+0.26 


-0.40 


+0.16 -0.39 


6 fragmentation 


±0 


71 


±0.56 


bjc semileptonic decays 


+0M 


-0.07 


±0.05 


Detector modeling: 








JES pt dependence 


±0 


25 


±0.19 


b response (h/e) 


+0.87 


-0.75 


+0.63 -1.43 


Trigger 


±0 


08 


+0.08 -0.13 


b tagging 






±0.24 


Method: 








Signal fraction 


+0.50 


-0.17 


±0.15 


QCD contamination 


±0 


67 


±0.29 


MC calibration 


±0 


17 


±0.48 


Total systematic uncertainty 


+1.5 


-1.4 


+1.2 -1.8 


Total uncertainty 


+5.2 


-7.5 


+4.3 -4.9 



TABLE X: Summary of uncertainties on the top quark mass. 
All values are quoted in GeV. 



Consequently, both calculations use a leading order 
parton distribution function (PDF): CTEQ5L [18]. 
To study the systematic uncertainty on mtop due to 
this choice, the variations provided with the next- 
to-leading-order PDF set CTEQ6M [21] are used, 
and the result obtained with each of these varia- 
tions is compared with the result using the default 
CTEQ6M parametrization. The difference between 
the results obtained with the CTEQ5L and MRST 
leading order PDF sets is taken as another un- 
certainty. Finally, the effect of a variation of Ug 
is evaluated. In all cases, a large ensemble has 
been composed of events simulated with CTEQ5L, 
and these have been reweighted such that distri- 
butions according to the desired PDF set are ob- 
tained. The individual systematic uncertainties are 
added in quadrature. The systematic uncertainty is 
dominated by that from the variation of CTEQ6M 
parameters. 

• b fragmentation: While the overall jet energy 
scale uncertainty is included in the statistical un- 
certainty from the fit, differences in the 6/light 
jet energy scale ratio between data and simulation 
may still affect the measurement. Possible effects 
from such differences are studied using simulated 
ti events with different fragmentation models for b 
jets. The default Bowler [22] scheme with r;, = 1.0 
is replaced with ri, — 0.69 or with Peterson [23] 
fragmentation with eb — 0.00191. Simulation stud- 
ies show that the variation of rt results in a change 
of the mean scaled energy (x_b) of 6 hadrons that is 
larger than the uncertainties reported in [24] , while 
the uncertainty on the shape of the xb distribu- 
tion is taken into account by the comparison of the 
Bowler and Peterson schemes. One large ensemble 
is built using events from each of the three simula- 
tions. The absolute values of the deviations of top 
mass results from the standard sample are added 
in quadrature and symmetrized. 

• b/c semileptonic decays: The reconstructed en- 
ergy of b jets containing a semileptonic bottom or 

charm decay is in general lower than that of jets 
containing only hadronic decays. This can only be 
taken into account for jets in which a soft muon is 
reconstructed. Thus, the fitted top quark mass still 
depends on the semileptonic b and c decay branch- 
ing ratios. They have been varied by reweighting 
events in one large ensemble of simulated events 
within the bounds given in [15]. 



B. Detector Modeling 

• JES Pt and |?7| dependence: The relative dif- 
ference between the jet energy scales in data and 

Monte Carlo is fitted with a global scale factor, 
and the corresponding uncertainty is included in 
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the quoted (stat. + JES) uncertainty. Any discrep- 
ancy between data and simulation other than a 
global scale difference may lead to an additional 
uncertainty on the top quark mass. To estimate 
this uncertainty, the energies of jets in the events 
of one large ensemble have been scaled by a factor 
of (1 + 0.02^0^^) where E-^^t is the default jet 
energy. The value 0.02 is suggested by studies of 
7+jets events. The top mass result from the mod- 
ified ensemble has been compared to the default 
number, and the symmetrized difference is taken 
as a systematic uncertainty. 

Similarly, to estimate the effect of a possible \ri\ 
dependence of the jet energy scale ratio between 
data and simulation, the jet energies have also been 
scaled by a factor (1 — 0.01|ry|) as suggested by 
7-|-jets events. No significant effect on the top 
quark mass has been observed and thus no addi- 
tional systematic uncertainty is assigned. 

• Relative 6/light jet energy scale: Variations 
of the h/e calorimeter response lead to differences 

in the 6/light jet energy scale ratio between data 
and simulation in addition to the variations of the h 
fragmentation function considered in Section VII A. 
This uncertainty has been evaluated by scaling the 
energies of h jets in one large ensemble and studying 
the effect on the top quark mass. 

• Trigger: The trigger efficiencies used in composing 

ensembles for the calibration of the measurement 
are varied by their uncertainties, and the uncertain- 
ties from all variations are summed in quadrature. 

• h tagging: The 6-tagging efficiencies are varied 
within the uncertainties as determined from the 
data, and the variations are propagated to the final 
result. 

Note that no systematic uncertainty is quoted due to 
multiple interactions/uranium noise as opposed to the 
Run I measurement. The effect is much smaller in Run II 
as a consequence of the reduced integration time in the 
calorimeter readout. It is moreover covered by the jet 
energy scale uncertainty, as the offset correction is com- 
puted seperately for data and Monte Carlo in Run II, 
accounting for effects arising from electronic noise and 
pileup. 

C. Method 

• Signal fraction: The normalization procedure 
of the background probability described in Sec- 
tion IV D is chosen such that the signal fraction 
/top as measured with the topological likelihood fit 
and given in Table II is reproduced. However, the 
signal fraction is slightly overestimated for low true 
signal fractions, which leads to a small bias in the 



resulting top mass. The signal fraction in ensem- 
ble tests used for the calibration is varied within 
the uncertainties determined from the topological 
likelihood fit, and the resulting variation of the top 
quark mass is taken as a systematic uncertainty. 

• QCD background: The VT+jets simulation is 
used to model the small QCD background in the 
selected event sample in the analysis. The system- 
atic uncertainty from this assumption is computed 
by selecting a dedicated QCD-enriched sample of 
events from data by inverting the lepton isolation 
cut in the event selection. The calibration of the 
method is repeated with ensembles formed where 
these events are used to model the QCD back- 
ground events whose fraction is given in Table II. 
The resulting change is assigned as a systematic 
uncertainty. 

• MC calibration: The calibration of the top mass 

measurement is varied according to the statisti- 
cal uncertainty of the calibration curves shown in 
Figs. 7 and 13. 



VIII. SUMMARY 

A measurement of the top quark mass using lep- 
ton-hjets tt events in 370 pb~^ of data collected with the 
DO detector at Run II of the Fermilab Tevatron Col- 
lider has been presented. The events are analysed with 
the Matrix Element method, which is designed to make 
maximal use of the kinematic information in the selected 
events. To avoid a large systematic uncertainty, an over- 
all scale factor JES for the energy of calorimeter jets 
is determined simultaneously with the top quark mass. 
This in-situ calibration of the jet energy scale helps re- 
duce the overall uncertainty on the top quark mass when 
combining with other measurements. 

The resulting top quark mass is 

mtop = 169.2t^:° (stat. + JES) t\l (syst.) GeV (42) 

for an analysis that uses only topological information, 
and 

mtop = 170.3l4:5 (stat. + JES) t\l (syst.) GeV (43) 

when 6-tagging information is included. The jet energy 
scale is JES = 1.048lg;Q4Q (stat.) in the topological anal- 
ysis and JES = 1.027lo;o32 (stat.) when b tagging is in- 
cluded, indicating consistency with the reference scale. 
The two results are consistent with each other. To obtain 
a value for the top quark mass in combination with other 
measurements, the second, more precise value should be 
used. 
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